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We compute the magnetization of graphene in a magnetic field, taking into account for generahty 
the possibihty of a mass gap. We concentrate on the physical regime where quantum oscillations 
are not observed due to the effect of the temperature or disorder and show that the magnetization 
exhibits non-linear behaviour as a function of the applied field, refiecting the strong non-analyticity 
of the two-dimensional effective action of Dirac electrons. The necessary values of the magnetic 
field to observe this non-linearity vary from a few Teslas for very clean suspended samples to 20-30 
Teslas for good samples on substrate. In the light of these calculations, we discuss the effects of 
disorder and interactions as well as the experimental conditions under which the predictions can be 
observed. 

PACS numbers: 75.70.Ak , 73.22.Pr 



I. INTRODUCTION 

The physics of graphene has attracted a huge 
amount of interest since its discoverj*i due to its 
unique physical properties as well as its potential 
technological importance. The electronic properties 
and its behavior in a strong magnetic field have 
been the focus of wide activity as recent reviews 
summarize^i^. 

The Dirac-like spectrum of clean graphene has 
been established, see e.g. Ref.lJ| and references 
therein. Moreover, the partition function of two- 
dimensional massless systems are known to exhibit 
strong non-analytic behaviour as a function of ex- 
ternal fields^i^. This kind of behaviour can lead, 
for example, to non-Ohmian conductivity due to 
Schwinger pair production at certain conditions^i^ 
or signatures of quantum criticality^. In the case of 
conventional metals, the magnetism is a result of two 
contributions, coming from the spin (Pauli contribu- 
tion) or the orbitals (Landau diamagnetism). In this 
work, we examine in detail the orbital magnetization 
of graphene in a magnetic field, including for gener- 
ality the possibility of a mass gap, and show the 
appearance of a non-linear dependence on the ap- 
plied field, as a consequence of non-quadratic field- 
dependence of the partition function. The Pauli 
magnetization of graphene is linear as a function of 
doping and much smaller than the orbital magne- 
tization especially for low carrier densities (chemi- 
cal potential close to 0)^. The possible nonlinear- 
ity due to the Pauli contribution, which is a known 
phenomenon, a consequence of the saturation of lo- 
calized moments at higher fields, is discussed in the 
part of this work where the experimental verification 
is proposed. 

The nonlinearity happens as a result of the failure 
of the linear response approximation, which is ap- 
plicable when there is a mass scale, larger than the 
applied perturbation. This scale controls the per- 
turbative expansion. For the case of graphene this 
could be e.g. a mass gap, the temperature or the 



impurity scattering rate. Since the magnetic energy 
scale in graphene is known to be exceptionally large 
due its linear dispersion, even at moderate magnetic 
fields the magnetic energy scale becomes dominant, 
and therefore it violates the basic assumptions of the 
linear response regime. 

The non-linearity survives even at relatively small 
magnetic field for sufficiently clean samples, there- 
fore it is not sufficient to compute only the zero 
field magnetic susceptibility. Typically, the observed 
non-linear effect of strong magnetic field is the mag- 
netic oscillations, but these oscillations average out 
at elevated temperatures or due to impurities pro- 
ducing the linear magnetization. Here we focus on 
the regime, where the magnetic oscillations aver- 
age out, but, nevertheless, a smooth non-linearity 
remains. In the following sections we discuss the ef- 
fects of temperature, mass gap, disorder, as well as 
interaction effects at a phenomenological level and 
propose an experiment which can demonstrate these 
predictions. 



II. CLEAN GRAPHENE WITH POSSIBLE 
MASS GAP. 

At zero temperature, the grand canonical poten- 
tial of clean graphene in a magnetic field when Lan- 
dau levels (LLs) are formed, at zero chemical poten- 
tial, with energy of the LLs En = s\gn{n) yj a\nB\ 
where a = 2hevp and n is an integer reads 



Qyac ^ gsgvMBl^'^^C 
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with a ~ 0.142 nm being the distance between the 
carbon atoms and the degeneracy factor C = ^^, 
gy = 2, gs = 2 accounts for spin and valley degen- 
eracy, C is the zeta function: '^^J^ ' = — C(— 1/2) ~ 
0.2079. 

The subleading terms are lattice corrections which 
are discussed elsewher o^°-^^ , the numerical fac- 
tor corresponds to nearest-neighbour tight-binding 



model calculation and describes orbital paramag- 
netic contributions due to higher energy levels hav- 
ing non-Dirac dispersion. The corrections as a result 
of the deviation of the dispersion from the Dirac- 
like appear either due to interaction effects - which 
we discuss below - or due to lattice effects when 
B ^ 100 Tesla and are not discussed further here. 

The leading term, though, is non-analytic in the 
magnetic field and naturally leads to divergent dia- 
magnetic susceptibility at the Dirac point. Such 
non-analyticities are typical not only in the effec- 
tive action of massless systems but at quantum crit- 
ical pointsi^rJ^ or even in corrections to Fermi liquid 
theorjii^. 

When temperature-averaged over many Landau 
levels (LLs), this non-analytic contribution cancels 
out, as first shown by McClure^. When the mag- 
netic field overcomes the energy scale set up by 
the temperature (or impurity scattering), this non- 
linearity can be observed. For generality and to com- 
pare energy scales, we take into account a mass gap 
A in the dispersion relation, which can be experi- 
mentally created e.g. due to A-B sublattice asym- 
metry caused by SiC substrate or by regular deposi- 
tion of impurities^SiSi. At zero magnetic field it was 
demonstrated that the pseudo-spin degree of free- 
dom (due to valleys) produces diamagnetic suscep- 
tibility which is^^"— : 



X(e) = -9sgv 



iC 



f2|A| 



(1) 



it is evident that the mass gap resolves the formal 
(5-function singularity in the susceptibility, but when 
the magnetic energy scale exceeds that of the gap, 
the non-analyticities of the free energy become again 
dominant. 

In a magnetic field B and in the presence of the 
mass gap, the spectrum at low energies becomes 



where we have regularized and subtracted the B ~ 
expression, C is a generalized Hurwitz ^-function. It 
also agrees with related formula in Ref. [22]. The 
above expression is exact in the limit of free elec- 
trons and has a strong-field expansion valid for 
5^ = ^^gr < I (when the gap is less or compara- 
ble to the distance between the first LLs), leading 
to the strong-field expansion of the magnetization 

(m = o, r = o): 

Af (A, 0, 0) = gsC^^B] (3C(-I/2) + S + ^^S^ 

ams^ _ m^s^ + 0(58)) (4) 
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for example, for 



VTbKt)' ^ 

A = O.IeV this forinula works starting from approx- 
imately 7.5 Tesla, indicating the scale from which 
the non-linearities occurs. 

At T = and as A — > the range of the values 
of B where linear magnetization occurs decreases to 
zero. Therefore, the nonlinearity of magnetization 
gets stronger with the reduction of the gap and tem- 
perature. The susceptibility, following from Eq.® 
reads: 
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At finite temperature and chemical potential, af- 
ter performing the same subtraction as for the reg- 
ularization of r2(A,/i = 0,T = 0), and additionally 
subtracting fJ-N(fi = 0, T = 0) which is independent 
of B we obtain^^: 



17(A,Ai,r) = f](A,0,0) 



(6) 
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E^^o = sign{nWa\nB\ + A^ ; Eo = (-1)"A (2) _^^ y ^^^ A ^ ^^^ f -E^+f^ 



where w = 0, 1 enumerates the two Dirac valleys. 
In the absence of the gap or when its size is smaller 
than the distance between the first LLs A < ■\/Q^p3| , 
the linear expression ([T]) is not applicable and, as we 
will show, the correct result leads to the non-linear 
magnetization. 

The regularized free energy reads (at chemical po- 
tential fi — 0, temperature T = 0) 



- [m -^ -Ai] 



At zero chemical potential the temperature plays 
a role similar to the gap. It is instructive to note 
that if we compare McClure's result for magnetic 
susceptibility at finite temperature: 
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with the one for the magnetization with gap but at 
zero temperature, we see that at 2kBT = A the 
linear magnetization is the same (with /x = 0), while 
the non-linear parts are different as is shown in Fig. 

m 

For non-zero chemical potential and temperature 
there are two regimes, the well known low-field 
regime and the high-field regime which sets in when 
El — \Ja\B\ > 2/i. When the chemical potential 
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FIG. 1. The magnetization at T=300 K, ^I = and at 
gap values A = and 52 meV is compared to the hn- 
ear behaviour. In addition, the magnetization at fixed 
chemical potential 100 meV and at fixed electron con- 
centration lO^^cm"^ at T=300 K is plotted. Low and 
high field asymptotic are shown. 



exceeds the temperature scale, the magnetization 
rapidly decreases as expected, but it starts to grow 
when the separation between the LLs becomes com- 
parable to the chemical potential, this behaviour is 
shown by the upper curve in FigH] At such fields the 
zeroth LL gives the leading constant paramagnetic 
contribution to the magnetization, but the same 
non-linear vacuum energy contribution remains: for 
V'^ > 2(^ -H A) > fcsT 

M(A, ^l, T) « g,C (// - A) + M(A, 0, 0) (7) 

this asymptotic regime is shown by the upper green 
dot-dashed line in Fig. [1] When the temperature 
decreases we get de Haas - van Alphen oscillations 
for non-zero chemical potential, as expected^^. The 
non-linearity we discuss can be alternatively inter- 
preted as being connected to this dHvA oscillations, 
as its remnant behaviour. 

When the number of particles is fixed, instead of 
the chemical potential (both situations are experi- 
mentally realizable in cases of graphene on substrate 
or suspended flakes) then ^ is expressed through the 
relation N = —^. At small temperature, the De 
Haas van Alphen oscillations are observed. At high 
fields the chemical potential inevitably tends to zero 
since almost all the electrons (or holes if /i < 0) can 
be hosted by the zero LL. This leads to the |i?|^''^ 
asymptotic behaviour of magnetization, correspond- 
ing to Eq. (I3|) , shown by the black solid curve in Fig. 

HI 

For completeness, we briefiy comment on the 
case of bilayer graphene where the zero-field 
susceptibility^SiSl diverges logarithmically with the 
Fermi energy e^ ^ and the divergence is cut by 
the greatest of trigonal warping scale ctrig and ep- 
When we incre ase th e magnetic field, the magnetic 
energy scale y^ap^l eventually becomes the great- 
est, thus leading to weak logarithmic non-linearity 



of the magnetization. The asymptotic form of 
the magnetic thermodynamic potential is^^: il — 

^^log(7i/(\MsT))S' where 71 « OAcV is 
th e inte rlayer hopping energy and the magnetic scale 
•\/ap3| is assumed to be larger than the trigonal 
warping energy etrig and ep, otherwise one replaces 
■\/ap3| — > ep . Then one gets 
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(8) 



At even larger magnetic fields B > 100 T, the mag- 
netic energy becomes larger than the interlayer hop- 
ping 71, thus effectively reducing the bilayer to two 
monolayers. Since the non-linearity of magnetiza- 
tion of graphene bilayer is significantly weaker than 
for the monolayer, we expect that the impurities 
would make this effect hard to observe. So, we con- 
centrate on the monolayer in what follows. The N- 
layered graphene was shown to have [iV/2] bilayer 
bands and N mod 2 monolayer bands^^. 



III. EFFECT OF IMPURITIES 

Besides the mass gap and temperature, the non- 
linearity of the inagnetization is influenced by impu- 
rities and interactions. We consider first the short- 
range scattering impurities and then the effect of 
charge inhomogeneities: the electron and hole pud- 
dles. 



A. Short-range impurities 

Consider for simplicity the short-range impuri- 
ties with momentum independent scattering. It 
is sufficient to adopt the self-consistent Born ap- 
proximation (SCBA)^i^. The treatment of the 
vacancy-type impurities^, or the commonly used 
phenomenological Lorentzian broadening leads to 
similar conclusions qualitatively, as we have checked. 

In SCBA^S, the self-energy reads 



S(e) 
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g{(^n) 



E(e) 



(9) 



where e„ is a full spectrum and g{e) is some smooth 

cut-off function. This is valid for small strength of 

2 
disorder W — ' j <^ 1. where Ui is the impurity 

concentration and Ui a measure of the strength of the 
random on-site impurity interaction. For discrete LL 
spectrum and for small W we obtain the solution 
iteratively. The density of states is 
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(10) 



For the strong-field regime the single-level 
approximation of the level width works well. 



as it has been checked numerically. Solving 
the Eq.® with a single energy level we obtain 
two roots and the relevant solution is: S(e) = 

i (e - e„ - sign(e - e,,) ^ {e - e,,)-^ - 2aW b) , 

therefore the density of states can be approximated 
by a semi-circle form 
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^/2aW\B\ - £2 0{2aW\B\ ~ e^) (11) 



where is a step-function and e means the deviation 
from the LL of the clean system. In the single-level 
approximation, we neglect the shift of the level cen- 
tre due to the real part of the self-energy. Similar 
form for the density of states is supported by more 
elaborate computation s, se e e.gi^. Note that the 
level width scales as \/|B|. The above approxima- 
tion can be applied for levels e^ up to n < gij, for 
higher levels their overlap would become essential. 
Since at medium or strong fields the higher levels 
are typically far from the Fermi surface, we may use 
the same level-broadening for all the levels, since this 
does not alter the result for the magnetization. 

Under the assumption that all the levels are 
broadened with the same profile p(e), the partition 
function can be computed by integrating Eq. Q with 
the broadened chemical potential: 



^ ^irnp / ^^ 



P(e-M) 



/p2(e')d< 



jn{A,e,T) (12) 



We note that this equation automatically inherits 
the regularization from il(A,e,T). For the magne- 
tization one gets a similar formula, but with an ex- 
tra contribution coming from the _B-derivative of the 
broadening profile. At high fields this contribution 
is paramagnetic due to broadening of the zero LL. 
As mentioned above, there is no significant depen- 
dence on the actual width of higher LLs. The effect 
of temperature may be alternatively taken into ac- 
count, by convolving the impurity-broadening with 
the derivative of the Fermi function with respect to 
energy /— . 

The typical line width of the high-quality 
samplo^f^ is estimated to be i5 w 3-\/i3(T) meV, 
so VF ss 0.003 and the Fermi energy is /i w 14meV. 
Such high-quality sample is close to the ideal case 
and exhibits strongly non-linear ^ -\/|i?| magnetiza- 
tion already at 1 Tesla, as illustrated in Fig. [5] 

From this subsection we conclude that low con- 
centrations of short-range impurities do not signif- 
icantly alter the magnetization. For larger concen- 
trations of impurities the above analysis is not ap- 
plicable. 



B. Charge inhomogeneities 

It was shown experimentally, that for many 
jraphene samples the charge inhomogeneities play 




B(T) 

FIG. 2. Non-linear magnetization with short-range im- 
purities for fixed fi or for fixed A'' . Green dashed curve 
takes into account the linear in B splitting of zero LL 
due to interactions, while the green dot-dashed curve is 
for ~ ^/B splitting of zero LL - see Section [TV] 



the dominant role^. The effect of charge inhomo- 
geneities can be modelled as long-range smooth vari- 
ation of carrier density A^ (electron and hole pud- 
dles), and the range may be assumed to be compa- 
rable to the size of crclotron orbits. A simple model 
proposed in Ref. [38| fits well the experiment. It 
assumes Gaussian variation of A'^ with standard de- 
viation of order 6N « 4 • lO^^cm"^. Since N ^ fi^, 
the variation SN leads to large 6fj. close to the tip 
of the Dirac cone and results in the large broaden- 
ing of the zeroth LL. This contrasts with the equal 
profile broadening, coming from the scattering on 
short-range impurities. Apart from variation of the 
charge density, the charged impurities cause a usual 
broadening of the levels. We use the simple model 
of constant Lorentzian broadening, independent of 
the magnetic field, as compared to short-ranged im- 
purities, where we got the ^/B dependence. This 
does not significantly change the result and is partly 
justified by computations in Ref. [35[, where the B 
dependence was shown to be more shallow than y/B 
due to the fact that the screening increases with B. 
To perform calculations, we broaden the levels 
with Lorentz profile, getting the density of states 
Pe{E) then integrate it to find N{Ef) and it's in- 
verse Ep{N) and then contract with Gaussian den- 
sity fluctuation profile: 



P{N,N,SN) 



1 



2tt6N 



cxp 



2dN^ 



(13) 



For example, the density of states is given by 



p{N) = / dNp,{EF{N))P{N,N,6N) (14) 

The consideration of graphene with fixed total num- 
ber of electrons is equivalent to the situation of 
graphene on substrate and the particle number be- 
ing proportional to gate voltage. For suspended 
graphene (or exfoliated flake) one may imagine a sit- 
uation of the fixed local chemical potential at spots 




M (eV) 

FIG. 3. Density of states for 16 T (solid) and 10 T 
(dotted) as a function of chemical potential of neutral 
spot with gaussian carrier density fluctuation SN « 4 ■ 
lO^^cm"^ and Lorentzian broadening of levels with half- 
width e = 15 meV. 

where it touches a contact. Assuming that these 
spots are far from the charged impurity (neutral 
spot), the local electron density at such spots would 
coincide with the average value N, so we may use 
the same as above function N{Ef) to convert p{N) 
to p{N{fi)). The resulting density of states (DOS) 
as a function of the chemical potential is plotted 
in the figure [3] for magnetic fields 10 and 16 Tesla, 
6N « 4 • lO^^cm"^ and Lorentzian broadening with 
half-width e = 15 meV. The same figure but plotted 
against the electron density N (or gate voltage) can 
be found in Ref. [3al. For high temperatures, the 
level broadening is dominated by temperature, and 
for low temperatures - by impurities. 

For low temperatures and Lorcntz impurity broad- 
ening e the Eq.([6|) for the grand canonical potential, 
combined with Ea.([T2|) gives at T = 0: 

n{A,n,e) = n{A,0,0)+ (15) 

+ g.,CB [F(p + A) + 2 ^^=1 ^(a^ + En)] + [m ^ -p] 

where En = \/an\B\ + A2, 

F{e) = -(arctan(max(A, e)/e) — arctan(max(— A, e)/e)) 

(16) 




FIG. 4. Non-linear magnetization with charged im- 
purities for fixed average electron density. Bottom to 
top: iV = 0; 5 • lO"; 10^^; 1.5 ■ lO^^cm"^, temperature 
r = 300 K (black solid) and T = (red dashed); Lorentz 
broadening e = 15 meV, density fluctuation dispersion is 
SN = 4- 10"cm"^ A plot for smaller SN = 4- 10^" cm"^ 
with N = lO^^cm"^ and T = (blue dotted) is shown 
for comparison. 

Inset: Dependence on the impurity strength. Red 
dashed curve is the same as the bottom one on the main 
plot: A = 0, T = 0, e = 15meV, <5A = 4 • lO^^cm"^; 
green dotted: Lorentz broadening increased to e = 
30meV; blue dot-dashed: density fluctuation increased 
to (5iV = 6 ■ 10"cm-^ 
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and A is a large cut-off for Lorentzian broadening. 

To see the effect of charge puddles on magnetiza- 
tion, we compute the magnetization as a function of 
electron density N and then convolve with Gaussian 
profile: 



5A = 4-lCcm" 
fi = 100 meV, T 



M{N, SN, T) 



2TrSN 



(N-N)-^ 



dNM{N,T)e ^(^^ 



(17) 

this prescription follows from summation over sepa- 
rate puddles, A^ denotes the average electron density. 
The results for fixed average carrier number N 
and for fixed chemical potential p of neutral spot 
are shown in Figs. [U [S] From these plots we see 
that charge disorder observed in experiments plays 
an important role to smear the magnetic oscillations 



FIG. 5. Non-linear magnetization charged impurities 
for fixed chemical potential of neutral spot. Bottom to 
top: /i = 0; 50; 100; 150 meV temperature T = 300 K 
(black solid) and T = (red dashed); Lorentz broad- 
ening is e = 15 meV, density fluctuation dispersion is 
* A plot for smaller SN = 4- 10^°cm-^ 
(blue dotted) is shown for compar- 
ison. 

Inset: Dependence on the impurity strength. Red 
dashed curve is the same as on the main plot: fi = 
150 meV, T = 0, e = 15meV, 5A = 4 ■ lO^cm"^; green 
dotted: Lorentz broadening increased to e = 30meV; 
blue dot-dashed: density fluctuation increased to SN = 



(the plot for lower charged disorder shows clear mag- 
netic oscillations), as does also the temperature, but 
the resulting magnetization is still non-linear. The 
non-linearity of magnetization gets weaker with in- 
crease of both types of disorder, but in a different 
way - see the insets to Figs. HI O 



IV. EFFECT OF INTERACTIONS 

It is worth noting that the most drastic effect of 
interactions is the splitting of the zero LL into two 
levels separated by a new gap 2A (different from the 
initial A), when ^ « and the zero LL is not com- 
pletely filled^i^'^. A grows with magnetic field. 
The precise form of its dependence on the inagnetic 
field is an open question and depends on the sample. 
Generally, there are the following kinds of energy 
gaps as a consequence of the magnetic field: (i) lin- 
ear in B dependence, coming from the Zeeman spin 
splitting Az = ifisB w 0.11S(T)meV and from 
the potential pseudospin splitting (from Kekule-type 
distortion of the lattice) Axekuie ~ 0.2_B(T) meV, 
and (ii) interaction contributions, scaling as vB, as 
explained e.g. in Ref.[3| and references therein. 

The experimental data in Ref.[37[ contain sig- 
nificant uncertainty, allowing for various fits. For 
samples with high mobility 17000cm^/(Vs) the fit 
by \/B looks reasonable: A « fci/2 ^/l^l, with 
ki/2 ~ 3 nicV/vT. For samples with lower mobility 
a fit linear in B is within error bars: A sa kunear \B\ 
with kiinear ~ 0.8 mcV/T at half-filling. This is 7 
times larger than the Zeeman splitting, indicating a 
different mechanism. 

The splitting of the zeroth LL leads to an ex- 
tra paramagnetic contribution, as compared to 
Eq.® with A = 0, /i = 0. The correc- 
tion to thermodynamic potential is easily com- 
puted by replacing log2 ^ ^ log fl + e^^/^'^^^M + 

^ log ( 1 + e+'^/f'^'^-^' ] . For the linear dependence 

- — - we have in the single- 



on B, defining x - 
particle approximation: 



fcflT 



SM ~ CgskBT ( 2 log cosh — h x tanh ■ 



This is shown by the dot-dashed curve on Fig. [21 
The paramagnetic nature of this contribution can 
be understood as due to the reduction in energy of 
the filled half of the zeroth LL by interactions. The 
correction is quadratic for small x and linear at large 

X. 

For \/B splitting, the correction of the magneti- 
zation is 



5M = CgskBT[ 2 log cosh 



fci/2 yi^l 

knT 



1 , X 

— .T tanh — 

2 2 



with X 

The different fittings of the splitting of the zeroth 
LL as a function of B, lead to different B-dependence 
of the magnetization, which, in turn, can provide 
one more method to distinguish between the main 
possibilities. We note that still there is no fully satis- 
factory theoretical explanation of the level-splitting 



effect^!, and its microscopic explanation is beyond 
the scope of this work. 

V. EXPERIMENTAL PROPOSAL AND 
DISCUSSION 

The measurements of non-linear magnetization 
can be used to study various properties of graphene 
sample: magnetization is sensitive to the number of 
carriers, mass gap and disorder as well as number of 
layers. In particular, one can extract the magnetic 
field dependence of the interaction-induced splitting 
of the zero LL. One possible way to measure the 
nonlinear magnetization is to measure the magneti- 
zation of a suspended graphene flake with scanning 
SQUID microscopy. Alternatively, one may go to 
higher magnetic fields (up to 50 Tesla) with larger 
amount of lower-quality graphene samples, e.g. the 
graphene laminate, used in Ref.[40]. At high fields 
the cyclotron orbits would fit better in the small- 
sized flakes and one can neglect the boundary effects. 
In real samples there will be a significant non-linear 
part of the magnetization, coming from localized im- 
purity spins (Pauli contribution), but this effect can 
be fitted at lower magnetic fields or fields parallel 
to the surface, and consistently subtracted^. More- 
over, at magnetic fields greater or of the order of 10 
Tesla and at low temperatures T < 4 K, it is ex- 
pected that all the localized moments will come to 
saturation, therefore the detected nonlinearity will 
be a consequence of the orbital contribution. 

In conclusion, the two-dimensional nature and the 
linear spectrum of graphene are the necessary con- 
ditions to observe non-linear magnetization at ac- 
cessible magnetic fields. The underlying reason is 
the breakdown of the linear response theory due to 
the fact that the magnetic energy is the dominant 
in the system. The linear dispersion relation gives 
relatively large distances between the LLs and the 
two-dimensionality leads to absence of k^ dispersion 
of LLs. We have found that even at room tempera- 
ture and with moderate concentration of impurities 
the non-linearity should be revealed at about 10- 
20 Tesla, while with very clean suspended samples 
at low temperatures a lower value of magnetic field 
is sufficient. There are two types of non-linear be- 
haviour that are present: near half-filling, the mag- 
netization scales as \/B at higher fields, due to non- 
analyticity of the effective action; then at higher val- 
ues of the chemical potential, the magnetization is 
small and linear at small fields, while it increases the 
slope after the first LL crosses the Fermi energy, and 
at even higher magnetic fields one can also observe 
the ^/B behaviour. 
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